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Abstract 

We model an infinitely long liquid bridge confined between two plates chemically pat- 
terned by stripes of same width and different contact angle, where the three-phase contact line 
runs, on average, perpendicular to the stripes. This allows us to study the corrugation of a 
contact line in the absence of pinning. We find that, if the spacing between the plates is large 
compared to the length scale of the surface patterning, the cosine of the macroscopic contact 
angle corresponds to an average of cosines of the intrinsic angles of the stripes, as predicted by 
the Cassie equation. If, however, the spacing becomes of order the length scale of the pattern 
there is a shaip crossover to a regime where the macroscopic contact angle varies between the 
intrinsic contact angle of each stripe, as predicted by the local Young equation. The results 
are obtained using two numerical methods. Lattice Boltzmann (a diffuse interface approach) 
and Surface Evolver (a sharp interface approach), thus giving a direct comparison of two pop- 
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ular numerical approaches to calculating drop shapes when applied to a non-trivial contact 
line problem. We find that the two methods give consistent results if we take into account a 
line tension in the free energy. In the lattice Boltzmann approach, the line tension arises from 
discretisation effects at the diffuse three phase contact line. 

Introduction 

The contact angle, between the tangent to a drop and the solid surface that supports it, is an im- 
portant concept in many different applications such as coatings, detergency, printing, adhesives 
and dentistry.- The contact angle provides information about how a liquid spreads on a surface 
in a given solid-liquid-gas system and it allows an estimation of the surface energy of solids.^ A 
drop of liquid on an ideal solid surface, neglecting gravity effects, will be a spherical cap with a 
circular contact line between the three phases and the same contact angle at all points around the 
contact line. However, topographic and chemical defects on the surface can lead to a contact line 
that is not circular, and a contact angle that varies along the contact line.- Traditionally such a 
variation, has usually been ignored as, for drops much larger than any surface feature, an effective, 
average contact angle is measured regardless of the observation direction. However, now that it is 
relatively easy to design well-defined micropatterned surfaces, and observe the behaviour of drops 
with dimensions of order the surface patterning,-i^i^i^i^i^ii^ variations in contact angle around the 
drop can be substantial and can be measured. 

We distinguish between two very distinct behaviours of a three-phase contact line on a pat- 
terned surface which can affect the uniqueness of the contact angle.— If the boundaries between 
regions of different wettability are perpendicular to the contact line as, for example, in Figure[2a), 
then the interface will adopt a corrugated shape to minimise its free energy and the contact an- 
gle will vary along the contact line. Such variation is known as contact angle multiplicity. This 
configuration corresponds to thermodynamic equilibrium and the final drop shape is independent 
of the initial conditions. 1 2, 13 ^14^15 however, the boundaries are parallel to the contact line, the 
drop can jump or be pinned leading to contact angle hysteresis. ■^^'^^^'^^^'^^^ Now the final drop 
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state depends on its dynamic history. An example of surface patterning where this behaviour will 
dominate is depicted in Figure [Hb). A drop placed at the centre of the axial pattern will spread 
until it is pinned by a (relatively) hydrophobic circle - which particular circle will be selected by 
the initial volume and energy of the drop, and will in turn determine the measured value of the 
contact angle. 



(a) (b) 




Figure 1 : A chemically-patterned substrate that will lead to a three-phase contact line that is (a) 
corrugated, but not pinned, (b) pinned, but not corrugated. 

In general, on real surfaces, both multiplicity and hysteresis in contact angle will be important, 
and the presence of one is typically accompanied by the other. However, there are experimental 
and theoretical works in the literature^>^'^2j2^ which have considered geometries aiming to sep- 
arate the two effects, and this is a helpful way to investigate a complicated problem. We follow 
this approach here, concentrating on modelling the equilibrium state of a liquid bridge confined 
between two chemically striped plates such that the contact lines run, on average, perpendicular to 
the stripes. We describe the crossover between the behaviour of the contact line and the contact an- 
gle when the spacing between the plates is large, or small, compared to the length scale of surface 
patterning. 

A second aim of our work is to compare two different numerical approaches, a diffuse interface 
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model, solved using a lattice Boltzmann code,— >^ and Surface Evolver,— an algorithm which 
assumes a sharp interface. We discuss the effects of the finite thickness of the interface and compare 
the efficiency and applicability of the two algorithms. 

In next section we describe the geometry of the model, outline the diffuse interface and Surface 
Evolver approaches and list the parameters used in the simulations. The results are then displayed 
and discussed. In particular, we explain how the surface patterning affects the macroscopic contact 
angle. Next, we summarise the paper and compare the two numerical methods. 

Methods 

Geometry 

An infinitely long liquid bridge was confined between two chemically-patterned walls as shown in 
Figure[2l These walls were smooth planes perpendicular to the z-axis, lying at z = 1 and z = H +\, 
and infinite in the y-direction. Both walls were patterned with stripes of equal width (A/2), lying 
parallel to the x-axis. Stripes on the two walls faced each other. The intrinsic contact angles of 
the stripes were taken to alternate between 0,- = 60° and 0,- = 30°. The width D of the liquid 
bridge, i.e. the bridge dimension parallel to the stripes, was fixed sufficiently large that there was 
no interaction between the two liquid-gas interfaces {D^ For the diffuse interface model, the 
interface thickness ^ provides an additional length scale. The effects of gravity are neglected in 
this paper. 

The interface lies, on average, parallel to the y-axis, with an oscillation because the fluid prefers 
to wet the stripes of lower contact angle, as shown in Figure [2tb). Our aim was to understand how 
the interface shape varies with the interface thickness, ^, the pattern period. A, and the spacing 
between the plates, H. 

We next summarise the two numerical approaches that were used to calculate the shape of the 
liquid-gas interface. These methods are two widely-used numerical approaches to calculating drop 
shapes. 
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(a) (b) 




Figure 2: (a) Liquid bridge between chemically-striped walls and (b) magnification of the liquid- 
gas interface at the three-phase contact line. 

Diffuse interface model 

The first numerical method is a mesoscale simulation approach where the equilibrium properties 
of the drop are modelled by a continuum free energy: 

^= / (v/^(n) + |(^«n)2)jy+ / V/,K)J5. (1) 
Jv 2. Js 

^fb{n) is a bulk free energy term which we take to be:— 

^fb{n) = PciVn + 1)2(V,? - 2V„ + 3 - 2/3Tvv) , (2) 

where V„ = {n — nc)/nc, Tw = {Tc — T)/Tc and n, ric, T, Tc and pc are the local density, critical 
density, local temperature, critical temperature and critical pressure of the fluid respectively. The 
parameter /3 is related to the density contrast between the liquid and gas phases. This choice of free 
energy leads to two coexisting bulk phases (liquid and gas) of density ^^.(l ± a/j8t^). The second 
term in Eq. ([T]) models the free energy associated with any interfaces in the system. The parameter 
K is related to the surface tension via y = {4-y/2Kpc{l5Xw)^/^nc) /3 and the interface thickness via 
= a/ Kn^. /AfiXwPc-— The final term in Eq. ([T]) describes the interactions between the fluid and 
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the solid surface. Following Cahn,— the surface energy density is taken to be V^5(n) = — ^n.v, 
where is the value of the fluid density at the surface. The strength of interaction, and hence the 
local intrinsic contact angle, 0/, is parameterized by the variable 0. In our simulations, chemically 
heterogeneous surfaces are simply modelled by setting the value of appropriately at every site of 
the solid surface lattice.— 

The dynamics of the drop is described by the continuity dB]) and the Navier-Stokes equations 

©: 

dtn + da{nua) =0, (3) 
dt{nua) + dp{nuaup) = -dpPap + vdp[n{dpUa + daUp + dapdyUy)] , (4) 

where u, P, and v are the local velocity, pressure tensor, and kinematic viscosity respectively. The 
thermodynamic properties of the system appear in the equations of motion through the pressure 
tensor P: mechanical equilibrium is equivalent to minimising the free energy. Eqs. [3] and |4] are 
solved using a Lattice Boltzmann algorithm which is described in detail in.— 1^^*22 

The liquid drop was initialised as a cuboid confined in the z-direction by the two chemically- 
striped surfaces, with periodic boundary conditions applied in the x and y directions. The simu- 
lation parameters which were used for all the numerical calculations were: K = 0.002, pc =1/8, 
He = 3.5, T = 0.4, Tc = 4/7, v = 0.1, and /3 =0.1, while those specific to a particular simulation 
are given at the appropriate place in the text. 

Surface Evolver 

Surface Evolver is a public domain software, developed by Kenneth Brakke,— which minimises 
the surface energy of a given volume of liquid within a prescribed geometry. The liquid-gas, 
liquid-solid and gas-solid interface energies, and hence implicitly any contact angles, are inputs 
to the model. If Surface Evolver is able to find the correct minimum as in,— >2^>^>^ it provides 
a useful alternative to diffuse interface models, both because it is computationally quicker, and 
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because in many cases a sharp interface represents the physically appropriate limit. 

To perform the Surface Evolver simulations the liquid drop was initialised as a cuboid confined 
in the z-direction by the two chemically-striped surfaces. Symmetry demands that the interface 
must meet the x — z plane at right- angles at the centre of each of the chemical stripes; we took 
advantage of this symmetry and imposed neutrally wetting walls at the centres of neighbouring 
stripes without altering the interface profiles. 

Measuring the contact angle and contact line corrugation 

Once the interface shape had been calculated, we recorded the interface profiles for different values 
of y. The macroscopic contact angle 0{y) was obtained by fitting a circle to the entire interface 
profile at y and measuring its angle of intersection with z = 1 . There will be a small correction 
because the contact line is not parallel to the j-axis, but this was found to be negligible. This 
definition of contact angle is similar to that typically used in experiments, e.g. in Axisymmetric 
Drop Shape Analysis (ADSA).— i^li^ 

The distortion of the contact line. Ax, was measured as the distance between the maximum and 
the minimum values of its x coordinate which occur, by symmetry, in the centre of the hydrophilic 
and hydrophobic stripes respectively. When the interface was diffuse, we defined its position as 
that where the density took the mean of its values in the liquid and gas phases. 

Results 

We aim to understand how the contact line corrugation, A)c, depends on the interface thickness, ^, 
the width of the stripes, A/2, and the height of the slab, H. To present the results we will scale all 
lengths to the spatial period A . 

We first consider the variation of the amplitude of the contact line distortion with the distance 
between the plates. There are two distinct regimes. For large bridge heights, the magnitude of the 
contact line distortion becomes independent of H/X. This occurs because the corrugations in the 
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interface decay with height over a healing length, of order X/2n,— small compared to the spacing 
between the plates. Hence the interface away from the surfaces is not corrugated. This is illustrated 
in Figure Oa), which shows cross sections across the liquid bridge at values of x corresponding to 
the centre of a hydrophobic stripe, the border between stripes and the centre of a hydrophilic stripe 
for H /X = 200/80. Figure [3tb) is a similar plot, but for plate separations H/X = 16/80. Now the 
decay length of the corrugation along z is larger than H and it is favourable for the corrugation to 
persist for all z. 

The variation of the macroscopic contact angles for different values of H /X, shown in Figure 
HI are a consequence of the behaviour described above. Figures I5a) and (b) were obtained using 
the diffuse interface model {^/X = 0.016) and Surface Evolver (without line tension) respectively. 
For large H/X, the macroscopic contact angle is independent of the position across the pattern, i.e. 
the contact angle multiplicity is mitigated. The value of the macroscopic contact angle is consistent 
with the Cassie angle, 6c, which corresponds to the arccosine of an average of the cosines of the 
intrinsic contact angles of the stripes.— The macroscopic contact angle obtained from the diffuse 
interface model deviates slightly from the Cassie angle (by ~ 2°). This deviation may be due to 
line tension effects and/or uncertainties in the simulation method. Drelich et. al.— pointed out 
that, when the contact line is contorted, there is a correction to the value of the effective contact 
angle due to the line tension. Simple estimates from our simulations show that the correction is of 
order 1°. Uncertainties in the lattice Boltzmann simulations arise from the discretisation errors in 
the implementation and measurement of the contact angle. This uncertainty is of order 2 — 3°. For 
the Surface Evolver data, the uncertainty comes from the measurement of the contact angle and is 
typically of order 1—2°. For small ///A, the macroscopic contact angle at the center of the stripes 
mirrors the local contact angle, while at the boundaries, it takes an intermediate value between the 
two local contact angles. 

Quantitative results showing the crossover between the two regimes described above are shown 
in Figure [5] which presents data for several different values of the reduced interface thickness, ^/X. 
The ^/X = results were obtained using Surface Evolver; the rest are from lattice Boltzmann 
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simulations. In each case we observe a similar dependence of corrugation on the distance between 
the plates. When the height of the channel is much larger than the healing length of the interface 
corrugation, the dimensionless deformation of the contact line saturates. This is as predicted by 
the classical theory of capillarity, because the contact line corrugation must be proportional to the 
characteristic length of the pattern.— As H/X is decreased there is a slight decrease in Ax/ A 
because the larger Laplace curvature between the plates inhibits corrugation. The typical value of 
Ax/?i obtained here is consistent with previous work by Hoorfar et. al— 

As H/X is decreased below 1, the distortion first decreases and then increases sharply. This 
increase occurs because the interface is in the regime where it remains corrugated for all values 
of z and smaller H reduces the excess interface free energy resulting from the corrugation, but not 
the wetting energy gained at the plates. In the inset of Figure |5] we show the Surface Evolver data 
close to the crossover region. 

To explore the variation of the corrugation with interface thickness more closely, Ac/ A is plot- 
ted against / A for three different values of /// A in Figure |6l For H/X = I and 1 the contact line 
distortion decreases slightly as the interface thickness becomes larger, reminiscent of the flattening- 
out effect that would result from including a line tension in the free energy. This suggests that the 
lattice Boltzmann model incorporates an effective line tension, due to discretisation effects, which 
increases slowly with increasing interface thickness. In the inset of Figured we plot Ax/X against 
^/(yA) obtained using Surface Evolver. ^ and 7 are defined as the line tension and the liquid-gas 
surface tension respectively. The Surface Evolver results show that Ax/X has a similar dependence 
on i^/y, as on ^ for the lattice Boltzmann simulations. This dependence of Ax/X on the line ten- 
sion is also consistent with previous studies by Neumann et. al.^ (and the references therein). 
For ///A = 0.2, however, the distortion increases slightly as the interface thickness becomes larger 
indicating that a diffuse interface favours corrugations along the solid surface. This behaviour is 
not reproduced in Surface Evolver when we take into account the effect of positive line tension. 
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— At the centre of the 60° stripe 

— At the boundary between stripes 

— At the centre of the 30° stripe 




Figure 3: Lattice Boltzmann results for the interface profile in the x — z plane at the centre of 
the 6i = 60° stripe, the centre of the 6i = 30° stripe and the boundary between stripes for (a) 
H/X= 200/80, (b) H/X = 16/80, and ^/A = 0.016. H, A, and ^ are respectively the height of 
the channel, the spatial period of the surface patterning, and the thickness of the interface. 



10 



Francisco J. Montes Ruiz-Cabello et al. 



Modelling the corrugation of the three- 



65 




0.5 1 1.5 2 2.5 3 



65 




2<i I I I I I I 

2.5 5 7.5 10 12.5 15 

Figure 4: Macroscopic contact angle at the centre of the 0,- = 60° stripe, the centre of the 0/ = 30° 
stripe and the boundary between stripes as a function of the height of the liquid bridge H. The 
interface thickness is set to (a) i^/A = 0.016 (diffuse interface model) and (b) i^/A = 0.0 (Surface 
Evolver). The error bars in (a) correspond to 3° and in (b) 2°. The lines are a guide to the eye. 
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Figure 5: Lattice Boltzmann results for the variation of the magnitude of the contact line corru- 
gation Ajc with the height of the liquid bridge H for different values of interface thickness . All 
lengths are scaled to the spatial period A and the lines are a guide to the eye. The inset shows the 
Surface Evolver data ((^/A=0.0) close to the crossover region. 
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Figure 6: Variation of the magnitude of the contact line corrugation Ax with the thickness of the 
liquid-gas interface for different values of the height of the liquid bridge H. All lengths are 
scaled to the spatial period A . The inset shows how the corrugation depends on the magnitude of 
the line tension C, in Surface Evolver. The line tension is scaled to the liquid-gas surface tension 7 
and the spatial period A. A is set to 0.8 and the lines are a guide to the eye. 
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Discussion 

We have presented numerical results for the behaviour of the interfaces bounding a liquid bridge 
confined between two plates patterned by stripes of differing contact angle for the particular case 
where the contact line runs, on average, perpendicular to the stripes. We were able to see clearly, 
using both a diffuse interface approach and Surface Evolver, a sharp crossover between a regime 
where the interface corrugations on the two surfaces are independent of each other, and decay 
moving away from the substrates, to a regime where the corrugation persists across the bridge. 
In the former case it is possible to define a unique macroscopic contact angle for the drop, as 
predicted by the Cassie equation. In the second the macroscopic contact angle varies between the 
intrinsic values on different areas of the surface, as predicted by the local Young equation. The 
crossover occurs for ///A of order unity. In both regimes, the interface diffuseness plays a relevant 
role through the interface thickness, which is related to the line tension. To compare the simulation 
parameters we have considered here to the physical variables, we note that the typical value of 
^/(yA ) in our paper is 10^^. Using 7 = 10^^ J/m^ and A = 100 nm— 1 /im, this corresponds to ^ 
= 10^^ — 10^^ J/m. This value of line tension is comparable to those reported in experiments.— 
Results from the diffuse interface algorithm approach those obtained using Surface Evolver 
(with line tension) in the limit that the interface thickness goes to zero. The advantages of Surface 
Evolver are that it is considerably quicker (typically by one or two orders of magnitude), that the 
line tension can be controlled, and that it immediately accesses the physical limit of an interface 
which is sharp on micron length scales. The advantages of diffuse interface models, on the other 
hand, is that they can model drop hydrodynamics, that they can follow changes in the topology of 
the liquid, such as drop break up, and that they can model pinning and depinning correctly beyond 
the quasistatic limit. It is pleasing that, for a problem such as this where both approaches should 
be applicable, the results for drop shapes are comparable. 
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